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Abstract. A model for positron binding to polar molecules is considered by 
combining the dipole potential outside the molecule with a strongly repulsive core 
of a given radius. Using existing experimental data on binding energies leads to 
unphysically small core radii for all of the molecules studied. This suggests that 
electron-positron correlations neglected in the simple model play a large role in 
determining the binding energy. We account for these by including polarization 
potential via perturbation theory and non-perturbatively. The perturbative model 
makes reliable predictions of binding energies for a range of polar organic molecules and 
hydrogen cyanide. The model also agrees with the linear dependence of the binding 
energies on the polarizability inferred from the experimental data [Danielson et al 
2009 J. Phys. B: At. Mol. Opt. Phys. 42 235203]. The effective core radii, however, 
remain unphysically small for most molecules. Treating molecular polarization non- 
perturbatively leads to physically meaningful core radii for all of the molecules studied 
and enables even more accurate predictions of binding energies to be made for nearly 
all of the molecules considered. 
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1. Introduction 

Positrons are a useful tool in many areas of science, such as condensed matter physics, 
surface science and medicine (see, e.g., nna). Despite this, there is still much about 
their interactions with ordinary matter that remains to be explored theoretically. In 
particular, the binding of positrons to matter has been a difficult subject to research 
[3]. On the part of theory, this is due to the strong electron-positron correlations 
which determine the binding energy and, in many cases, ensure the very existence 
of bound states. On the experimental side, positron binding to atoms has not been 
verihed experimentally, largely due the difficulty in obtaining the relevant species in 
the gas phase. On the other hand, for polyatomic molecules a wealth of information is 
now available thanks to the special role that vibrational Feshbach resonances play in 
positron-molecule annihilation [1]. 
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Before a positron annihilates with an electron in a molecule, it usually forms a 
quasibound state with the molecule by transferring its excess energy into vibrations of 
a single mode with near-resonant energy. This leads to pronounced resonances observed 
in the positron-energy dependence of the annihilation rate HE]. Using the relation 

€^=UJ^- Cb, (1) 

where is the energy of the resonance due to vibrational mode u with energy 
experimentalists have now been able to measure values of the positron binding energy 
eb for over sixty molecules BIZ]. These measurements led to the construction of a 
phenomenological parametric ht of eb in terms of the the dipole polarizability a and 
permanent dipole moment p of the molecule: 

eb = 12.4(a + 1.6/i-5.6), (2) 

where eb is in milli-electron volts, a is in cubic angstroms and p is in debyes (D) [8]. 
An interesting feature of ([2]) is that the dependences of eb on p and a are both linear. 
Although a general increase of eb with p and a is to be expected (since both contribute 
to the positron-molecule attraction), there is no obvious reason why the dependences 
should be linear. In fact, measurements for some molecules with large dipole moments, 
such as acetone and acetonitrile, yielded binding energies more than double the values 
predicted by equation ([2]) [6]. 

Despite the wealth of experimental data on positron-molecule binding energies, 
theoretical developments are somewhat behind. There are few calculations of positron 
binding to nonpolar or weakly polar molecules. The zero-range potential model [9l[in] 
captured the qualitative features of the binding for alkanes and correctly predicted the 
emergence of the second bound state [9]. There were also predictions of positron binding 
to the hydrogen molecule in the excited A state m, and conhguration interaction 
calculations for carbon-containing triatomics (CO 2 , CS 2 , CSe 2 and weakly polar COS, 
COSe and CSSe) [121 IS]. The latter papers reported binding by the two heaviest 
species in the vibrational ground state, by CS 2 in the lowest vibrationally excited states, 
and by other molecules at higher vibrational excitations or upon bond deformations. 
In contrast, there is a large number of quantum chemistry calculations of positron 
binding with strongly polar polyatomic molecules with dipole moments > 3D. For such 
molecules binding is achieved even at the lowest, static-held (e.g., Hartree-Fock) level of 
theory. The static-held binding energies, however, are usually quite small, and the ehect 
of correlations (e.g., polarization of the molecule by the positron) increases the binding 
energy dramatically (see, e.g., m Eg nEi ini E8]). Recent conhguration interaction 
calculations for nitriles, acetaldehyde, and acetone usi izniiii] in fact give binding 
energies within 25-50% of the experiment, which is quite good, given the complexity of 
the system. 

The purpose of this article is to present a simple model for positron binding to 
polar molecules. For many molecules of interest the dipole moment is dominated by 
a single bond (e.g., CN or CO), located at one end of the molecule, with the negative 
charge on the terminal atom. Given the positron repulsion from the atomic nuclei, we 
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model the molecular potential as a point dipole surrounded by an impenetrable sphere. 
Of course, the true size of the molecular dipole is finite, of the order of interatomic 
distances. However, for weakly bound positron states, the wave function of the positron 
is very diffuse. Its spatial extent is much greater than the physical size of the dipole, 
which justihes the applicability of the point-dipole model to weakly bound positron 
states. This is illustrated by figure [H which shows the density for the positron bound 
in the dipole held of the acetonitrile molecule (CH 3 CN). The hgure also shows that the 
positron is localized in the negative-energy well of the dipole potential and is largely 
“unaware” of the true geometry of the molecule. This justihes the hard-sphere model 
for the short-range positron repulsion. Quantum chemistry calculations of the positron 
density in polar molecules support the picture of a dihuse positronic cloud localized oh 
the negatively charged end of the molecular dipole [Ml [171 US EQ]. 


/?(au) 



Figure 1. The density of the positron bound in the field of a point dipole with the 
dipole moment y = 3.93 D of the acetonitrile molecule and repulsive core of the radius 
ro = 1.175 au, with the binding energy of Cb = 27meV. 

Note that a recent paper [22] combined a hard-sphere repulsive core with the 
polarization potential to model positron binding to atoms and nonpolar molecules. 
While the two models bear some similarity, the physics of positron binding to neutral 
atoms and nonpolar species is very diherent from that of binding to strongly polar 
molecules explored in this work. In the former case, for atoms the positron does form 
a spherical cloud, but for molecules the shape of the positron wave function largely 
repeats that of the molecule [lO], and the spherical repulsive core approximation is hard 
to justify. In contrast, for bound states with polar species, the positron resides in the 
dipole-generated well to one side of the molecule (see hgure H]), and the repulsive core 
model looks more appropriate. 

The main features of binding by the dipole potential are outlined in section [2j 
Applying the model to polar molecules for which the positron binding energies are 
known from experiment (section (Sj) shows that electron-positron correlations have a 
large effect on binding. We include these in the form of the polarization potential, hrst 
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via perturbation theory (section 0]) and then non-perturbatively (section 0]). While this 
model may appear to be rather crude, it captures the main physical aspects of the 
problem. Owing to its simplicity, it provides a deeper understanding of some of the key 
features that have been observed in experiment, including the empirical scaling ([2]). The 
usefulness of such models as a means of obtaining an explanation, and complementary 
to heavy numerical computations, was argued well by Ostrovsky, who introduced the 
notion of complementarity between calculation and explanation |23l |24] . 


2. Theory 


We model the molecule as an impenetrable sphere of radius ro, with a point dipole of 
dipole moment /r fixed at its centre (the origin). The positron experiences point-dipole 
potential in the region outside the sphere. Using spherical polar coordinates (r, 6, f) 
and choosing the polar (z-) axis along /x, we have 


Ki(r) 


oo for r < ro, 

pr~^ cosd for r > ro. 


( 3 ) 


where d is the polar angle, and we work in atomic units (an). 

Although ([3]) is a non-central potential, the Schrodinger equation. 


+ Vi(r) 


f’(r) = Eif(r), 


( 4 ) 


for the positron wave function if(r) and energy E can be solved in the region r > tq 
using separation of variables. Inserting the ansatz "^(r) = i?(r)<f)(0, 0) into (jlj) yields 
separate radial and angular equations: 


IT +1 

r^ dr \ dr / 

1 d f . 

Mm V 


, 1 

H—+ (A — 2/icos0)$ 

sm^ 9 


0 , 


(5 a) 
(56) 


where A is a separation constant. If /i = 0 then (1^ becomes {Lf — A)<h = 0, where 
Lf is the squared angular momentum operator. This is just the eigenvalue equation for 
the Lf operator; the possible values of A are l{l -|- 1), where I (the azimuthal quantum 
number) is a non-negative integer, and the eigenfunctions are the spherical harmonics 
Yim{9,(j)), where m (the magnetic quantum number) is the eigenvalue of L^, and m is 
an integer, \m\ < 1. 

For p 7 ^ 0, Hs no longer a good quantum number since does not commute with 
the Hamiltonian. However, does commute with the Hamiltonian, and thus m is still 
a good quantum number. We must solve the angular equation fl^ to find the new 
values of A. With this information we will be able to solve the radial equation and use 
it to investigate the dependence of the binding energy eb = \E\ on p and tq. 
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Since m is a good quantum number, there will be a distinct set of eigenfunctions 0) 

of the angular equation (1^ for each value of m. We expand the unknown functions in 
the basis of spherical harmonics, i.e.. 


^ (m = 0, ±1, ±2,...), (6) 

V=\m\ 

where the Cum are unknown numbers. Substituting this expression into ([56]), multiplying 
across by f), where I is a non-negative integer such that I > \m\, integrating over 

0 and 9 and using properties of spherical harmonics (see, e.g., IZ51) yields 


^ ^ Bll'mCl'ra — ^mClmi 


l'=\m\ 


(7) 


where 

Bii'm = l{l + 1 )^«' + 2 /i(-l)"*\/( 2 / 1 )( 2 /' 1 ) r V ( 8 ) 

\0 0 0/ \0 —m m J 

and the arrays in parentheses are 3j symbols. The eigenvalue A has been renamed Am 
since it will have different sets of values depending on m. Equations (Cj) are a set of 
matrix eigenvalue equations for the semi-inhnite, symmetric, tridiagonal matrix Buim, 
whose rows and columns are enumerated by I and I', respectively. Each of these matrices 
has a countably inhnite set of eigenvalues, so we rename Am as X^s, where s = 1, 2, 3,... 
enumerates the different eigenvalues for each m, and the eigenvalues are arranged so that 
Ami < Am 2 < Ams < .... Symmetry properties of the 3j symbols can easily be used to 
show that Biii _m = i?/z'm, and so o need only be solved for m > 0. 

We seek bound states of the positron. From the form of the radial equation (15^ 
it can be shown that for Ams < there will be an inhnite number of bound states, 
while for Ams > there will be none [26] (see section . Given a certain value of 
/i and of m, by truncating the inhnite matrix to a hnite size (where the hnal row 
and column are denoted by / = /' = /max), we can hnd numerical approximations for the 
hrst /max — |m| -|- 1 values of Ams- Table [Hand hgure[2]show how Ami, Am 2 and Ams vary 
with p for m = 0 and m = ± 1 . Values of /max shown in the last column of table [T] are 
chosen so that the eigenvalues are correct to at least six decimal places. 

Considering the eigenvalues Ams as functions of p, for each combination of m and s 
there is a critical dipole moment pcrit for which Ams = —j- Some of these critical dipole 
moments are shown in table [21 they agree with the values obtained by Fermi and Teller 
[27] and Crawford [28] . The condition p > p^rit guarantees binding by either a point-like 
or hnite dipole. The smallest critical dipole is pcrit = 1.625 D for m = 0, s = 1. Since 
typical molecules have dipole moments not exceeding 11D, i.e., up to 4.3an, the only 
possible bound states are those corresponding to m = 0, s = 1 and m = ±1, s = 1. The 
critical value of p needed to sustain any other bound state is simply too high. 

It must be mentioned that the above considerations apply to binding by the 
static dipole, i.e., assuming that the molecules cannot rotate. When rotations are 
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Table 1. Values of Xms for m = 0,±l,s = l,2,3 across a range of values of the 
dipole moment g. The values of /max needed for stability to six decimal places are also 
shown. 


|m| 

g (au) 

^ml 

Am2 


^max 

0 

0 

0.000000 

2.000000 

6.000000 

2 

0 

2 

- 1.704857 

2.602337 

6.412828 

7 

0 

4 

- 4.519910 

2.263955 

7.444429 

8 

0 

6 

- 7.616374 

1.031162 

8.141444 

9 

0 

8 

-10.856049 

-0.662826 

8.138189 

10 

0 

10 

-14.186766 

-2.640671 

7.597027 

10 

1 

0 

2.000 000 

6.000000 

12.000000 

3 

1 

3 

0.489 539 

6.153928 

12.285114 

8 

1 

6 

- 2.675 243 

5.535499 

12.865402 

9 

1 

9 

- 6.481474 

3.913801 

13.078672 

10 

1 

12 

-10.628924 

1.630021 

12.620 576 

11 

1 

15 

-14.995 686 

-1.093056 

11.558 712 

12 



Dipole moment p (au) 



Figure 2. Dependence of the eigenvalues Ami (purple circles), Am 2 (red squares) 
and Am 3 (blue triangles) for (a) m = 0 and (6) m = ±1, on the dipole moment g. 
Intersections with the dashed lines (A = —1/4) give critical dipole moments gait- 


included, the values of pcrit required for the dipole binding to occur are 10-30% greater 
|29] . This gap depends on the moment of inertia of the molecule and increases with 
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Table 2. Critical dipole moments required for various bound states of the positron, 
along with the values of ^max required for stability to 10“® au. 


|m| 

5 

/^crit (S-u) 

/^crit (D) 

^max 

0 

1 

0.639 315 

1.625 

4 

0 

2 

7.546 956 

19.182 

10 

1 

1 

3.791968 

9.634 

8 

1 

2 

14.112115 

35.869 

12 

2 

1 

9.529 027 

24.220 

10 


the molecular angular momentum, being smallest for large, slowly rotating molecules. 
Another consideration important for ab initio quantum-chemistry calculatons of binding 
is that for the values of p only slightly exceeding /icrit, the binding energy is very sensitive 
to the actual value of the dipole moment (see figure [3] in section 12311 . The actual value 
of the dipole moment depends on the approximation used (e.g., Hartree-Fock), and can 
be a significant source of error PI. However, both the effect of molecular rotations and 
the sensitivity to the value of p are offset by the large contribution of electron-positron 
correlations to the positron binding (see sections lH and [5l) . 


2.2. Radial equation 


Under the substitution 


R{r) = 


Z{kr) 

\/kr 


(9) 


where k = y/^E, the radial equation fl53) yields the following differential equation for 
Z{kr): 


,, ,2 d^Z , dZ 
[kr) ,,, -b kr- 


d{kr) 


d{kr] 


+ 


{krf - -b - 


= 0 . 


( 10 ) 


This is just Bessel’s differential equation. Since for bound states we have E < 0, i.e., 
Z = —1£'|, it is best to express the general solution in terms of modified Bessel functions: 


Zms{kr) AmsKi0rnsi.^'^) T EmsRlSrnsi^'^) y 


( 11 ) 


where the subscripts m and s have been added to Z because there is a distinct function 
for each combination of m and s, and B^s are arbitrary constants, k = -^2|i7|, and 

1 1/2 


f^ms 


ms + 4 


and Ar 


< 


•1 is assumed. Equation ([9]) then gives 


RmsiZ^ 


K; 


iff 

iPms \ / I D 

“ “T 


h 


i/3n 


Kr 


Kr 


( 12 ) 


(13) 


For the bound-state wave function to be normalizable we must require Rmsiz) —?■ 0 
as r —)■ cx). It can be seen from the asymptotic forms of the modified Bessel functions 
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(see, e.g., m) that /y(x)/y^ —)> oo as X —)■ oo, while Ky{x)/^ —)■ 0 as x —)■ oo, 
assuming that x is real. We therefore require B^s = 0 for every m and s, and so 




y/KT 


(14) 


with A^s a constant of normalization. 

Since nr and fdms are real and positive, the function (also known as the 

Macdonald function) is also real, which can be seen, e.g., from the integral representation 

1321 , 

poo 

= / exp(-fi;r cosht) cos(/9mst) dt, (15) 

Jo 

and thus Rmsi'f') is real (for a real Ams)- The function Rms{f') has inhnitely many 
positive roots, with an accumulation point at r = 0. 

The second boundary condition to be applied to Rms{r) is due to the impenetrable 
sphere at r = tq, which means that we must have Rmsijo) = 0, i.e.. 




= 0 . 


(16) 


The positive roots of the function therefore the allowed values of nrQ. Since 

these roots form an inhnite, discrete set, we shall name them (msn, where n = 1, 2, 3,..., 
and Cmsi > Cms 2 > Cms 3 > ■ ■ ■■ These roots can be found numerically. For any particular 
molecule, tq is a constant, and so the permissible values of k (which we now rename 
Kmsn, and likewise with E) are Kmsn = Cmsn/ro, i.e., 

.2 \ 2 






Cn 


ro 


( 17 ) 


2.3. Dependence of binding energy on ro and p 


For a given value of the dipole moment, the largest negative value of Xms is for m = 0, 
s = 1, with the critical dipole moment pcrit = 0.6393au = 1.625D [28]. The 
corresponding ground-state binding energy is 

2 


= - 
2 


1 ( Con 


ro 


( 18 ) 


where Con is the largest root of whose index /3oi is determined by the eigenvalue 

Aoi of the angular equation, see flT^ . The dependence of the binding energy fITHD on tq 
is simple. Figure |3] shows the dependence of Cb on the magnitude of the dipole moment 
p for the hxed repulsive core radius ro = lau (i.e., Cb = |Coii)- 

For /i —)■ /icrit, the binding energy rapidly tends to zero. This limit corresponds to 
Aoi —)■ —I and (3oi —t 0. In the limit of small /3, the roots of the Macdonald function 
Kijs^z) have the following asymptotic behaviour [33] : 
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Figure 3. Dependence of the ground-state binding energy Cb on /r for ro = 1. Solid 
red line with circles shows the numerical results, and the blue dashed line is the fit 
([231) with A = 65 998.7meV, B = 12 . 27050 ^/ 2 ^ and C = 0 . 400612 D- 1 / 2 . Note that 
the two curves are indistinguishable on the scale of the graph. 


where 7 ~ 0.577 is Euler’s constant. The largest root that we are interested in (n = 1) 
is then given by 


Co„^2e-'exp(-^). 

For dipole moments close to the critical value, we have from equation ffT2ll . 

1/2 


dA 


01 


d/i 


(/t /^crit) 


M—Merit 

and 


/3oi - 

Combining equations lli^p . li: 2 U|) and 11 : 21 |) . gives 
€b = /lexp , 


( 20 ) 


( 21 ) 


( 22 ) 


where A and B are constants (see also [3H |35], from which a similar result can 
be derived). Motivated by this scaling, we constructed an approximate analytical 
expression for the binding energy as a function of p in the following form: 

Cb = x4exp [-B{p - + C{p - • (23) 


Here the second term in the exponent represents a correction to the leading term 
(| 22 |) . It accounts for the next order corrections in both fl 20 ll and (l 2 TD . and extends the 
applicability of fIMll way beyond the range of near-critical p. Regarding the constants 
A, B and C as htting parameters, an excellent £t of the numerical data over the whole 
range covered by figure [3] is obtained using A = 65 998.7meV, B = 12.2705 and 
C = 0.400 612and the dipole moment p in debye (D). 
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Table 3. Values of rg obtained for a selection of molecules by fitting known binding 
energies from [3IM1 [3H] to equation (|T5)) . 


Molecule 

m(d) 

Cb (meV) 

ro (au) 

Aldehydes 




Acetaldehyde (C2H4O) 

2.75 

88 

1.03 X 10-1 

Propanal (CgHeO) 

2.52 

118 

4.28 X 10-2 

Butanal (C4H8O) 

2.72 

142 

7.44 X 10-2 

Ketones 




Acetone (CgHeO) 

2.88 

174 

1.00 X 10-1 

2-butanone (C4H8O) 

2.78 

194 

7.45 X 10-2 

Cyclopentanone (C5H8O) 

3.30 

230 

1.90 X 10-1 

Formates 




Methyl formate (C2H4O2) 

1.77 

65 

4.04 X 10-® 

Ethyl formate (C3H6O2) 

1.98 

103 

9.76 X 10-3 

Propyl formate (C4H8O2) 

1.89 

126 

1.77 X 10-4 

Acetates 




Methyl acetate (C3H6O2) 

1.72 

122 

7.38 X 10-8 

Ethyl acetate (C4H8O2) 

1.78 

160 

4.32 X 10-8 

Nitriles 




Acetonitrile (C2H3N) 

3.93 

180 

4.54 X 10-1 

Propionitrile (C3H5N) 

4.05 

245 

4.37 X 10-1 

2 -methylpropionitrile (C4H7N) 

4.29 

274 

5.04 X 10-1 

Methyl halides 




Methyl fluoride (CH3P) 

1.86 

0.3 

1.72 X 10-3 

Methyl chloride (CH3CI) 

1.90 

25 

1.68 X 10-3 

Methyl bromide (CHgBr) 

1.82 

40 

1.67 X 10-6 


3. Positron binding by the dipole potential 

Using experimental data on positron binding energies from [3136] and dipole moments 
from [ 33 , we fitted the energies to equation flT8|) by adjusting the values of tq for fourteen 
polar organic molecules and three polar inorganic molecule^. No direct experimental 
binding energy is available for methyl fluoride (CH 3 F), and the value obtained by htting 
theoretical annihilation rate to experiment has been used instead | 3 B] • If the static dipole 
potential provided the dominant contribution to the binding, then it could be expected 
that for molecules with a single dipolar bond the values of r^ would be approximately 
half the length of the molecular dipole (~1 an), with values signihcantly larger or smaller 
than this considered as unphysical. For molecules with several dipolar bonds (e.g., the 
formates and acetates) we expected a larger value of tq than in the case of a single dipole 
bond. 

The results are shown in table O Clearly, all of the radii obtained are unphysically 
small, particularly for the most weakly polar molecules. The largest value, r^ = 0.5 an, 

f The methyl halide molecules are quite distinct from the other molecules studied. Each of them 
contains a different halogen atom. It is this, rather than the size of the molecule, that affects their 
dipole polarizability. 
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is for the most strongly polar molecule studied: 2-methylpropionitrile, but even this is 
less than a quarter of the C=N bond length. 

It can be seen from table [3] that despite molecules of the same type (aldehyde, 
ketone, etc.) having similar dipole moments, there can be signihcant variations in the 
binding energies. For example, consider the molecules acetaldehyde and butanal. Their 
dipole moments are very close (2.75D and 2.72D, respectively), and the dipole in both 
molecules is due to a C=0 bond. Yet there is a large difference in the binding energies: 
the binding energy for butanal (142 meV) is more than 1.5 times that for acetaldehyde 
(88meV). Peculiarly, acetaldehyde is actually slightly more polar than butanal, yet has 
the lower binding energy. A similar situation also occurs with acetone and 2-butanone, 
and with ethyl formate and propyl formate. These observations cannot be explained by 
our model as it stands. 

The fact that the values of ro obtained for all of the molecules studied are 
unphysically small implies that lepton correlations (in particular, due to polarization 
of the molecule by the positron) play an important role in enhancing the binding 
energy, even for strongly polar molecules. Going back to the example of acetaldehyde vs 
butanal, acetaldehyde has a polarizability of 4.6 A^, while the polarizability of butanal is 
a signihcantly greater value of 8.2 [S^. This explains the larger binding energy of the 

latter molecule. In the following two sections we investigate the effect of the molecular 
polarization on positron binding, and show that its inclusion is critical for obtaining a 
correct physical picture of positron binding to polar molecules. 

4. Perturbative correction due to molecular polarization 
4-1. Core radii for perturbative inclusion of polarization 

Since the dipole potential for p > pcrit is sufficient to create a “zeroth-order” bound 
state, we hrst estimate the effect of molecular polarization using perturbation theory. 
The values of the radius ro can then be chosen by htting the total binding energy (i.e., 
due to the dipole force and polarization) to the experimental values, expecting that this 
should lead to more realistic values of vq. 

The extra contribution to the positron potential energy (in the region r > tq) due 
to molecular polarization can be approximated by the polarization potential 

^Poi(r) = (24) 

where a is the molecular dipole polarizability. Using perturbation theory, the hrst-order 
correction to the original dipole binding energy (IT^ . which we now label is 

= /^1^011 (r)|M^r, (25) 

where fimsnij) = Rmsn{r)^rns{6i<i>)- Assumiug that the radial and angular parts of the 
wave function are separately normalized to unity, this becomes 

poo 

=n \Roii{r)\‘^r~‘^ dr. 

^ J ro 


( 26 ) 
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Table 4. Fitted values of rg with the inclusion of polarization via perturbation theory. 
Also shown are the corresponding values of and , the expectation values of the 
potential energy due to the permanent dipole, and the predicted and experimental 
values of eb- 



h 

(D) 

a 

ro 

(an) 

4°^ 


l(Ki(r))| 

(meV) 

Cb (meV) 

Molecule 

(A3) 

(meV) 

b 

(meV) 

Pred. 

Exp. 

Aldehydes 

Acetaldehyde 

2.75 

4.6 

0.60 

3 

80 

38 

83 

88 

Propanal 

2.52 

6.5 

0.42 

1 

118 

23 

119 

118 

But anal 

2.72 

8.2 

0.58 

2 

140 

35 

142 

142 

Ketones 

Acetone 

2.88 

6.4 

0.63 

4 

168 

58 

172 

174 

2-butanone 

2.78 

8.1 

0.58 

3 

187 

46 

190 

194 

Cyclop ent anone 

3.30 

9.0 

0.92 

10 

220 

98 

230 

230 

Formates 

Methyl formate 

1.77 

5.1 

0.004 

~10-3 

69 

~io-2 

69 

65 

Ethyl formate 

1.98 

6.9 

0.066 

~10-2 

109 

2 

109 

103 

Propyl formate 

1.89 

8.8 

0.0305 

~10-3 

126 

~10-3 

126 

126 

Acetates 

Methyl acetate 

1.72 

6.9 

0.0006 

~10-6 

116 

~io-^ 

116 

122 

Ethyl acetate 

1.78 

8.6 

0.0048 

o 

1 

156 

~io-2 

156 

160 

Nitriles 

Acetonitrile 

3.93 

4.4 

1.175 

27 

155 

202 

182 

180 

Propionitrile 

4.05 

6.3 

1.24 

31 

218 

218 

249 

245 

2-methylpropionitrile 

4.29 

8.1 

1.40 

35 

244 

235 

279 

274 

Methyl halides 

Methyl fluoride 

1.86 

2.4 

0.042 

~io-^ 

0.33 

~10-2 

0.33 

0.3 

Methyl chloride 

1.90 

4.4 

0.026 

~10-3 

23 

~10-3 

23 

25 

Methyl bromide 

1.82 

5.6 

0.0085 

~10-3 

42 

~10-3 

42 

40 


For each of the molecules studied we made an initial estimate for the value of uq 
and adjusted it until the new binding energy Cb = 6^°^ + was within 10% of the 
experimental value (with both and being functions of tq). The results are shown 
in table m Polarizabilities are taken from the CRC Handbook [3^, with the exceptions 
of propyl formate and cyclopentanone, for which the polarizabilities have been estimated 
by Danielson ef a/ [7]. 

All of these new values of rg are signihcantly larger than the original values. The 
three nitriles are the most strongly polar molecules studied, and they now have very 
realistic values of tq. The dipole in these nitriles is due to the C=N bond. The length 
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of this bond is 2.19an [3S], half of which is approximately 1.1 an. The values of uq are 
only slightly greater than this. 

The ketones — acetone, 2-butanone and cyclopentanone — are the second most 
polar group. The polarity of these molecules is due to a C=0 bond, the length of which 
is 2.26 an [39] (half of this is 1.13 an). Their values of ro are not as close to this estimate 
as those for the nitriles. For the most polar molecule in the group, cyclopentanone, tq 
is within 19% of half the bond length. The values of rg for acetone and 2-butanone 
are, however, significantly smaller. The picture is similar for the aldehydes, which show 
vq ~ 0.5 an. The three other groups (formates, acetates and methyl halides) have dipole 
moments /i < 2D, only slightly exceeding the critical dipole moment /icrit = 1.625D. 
They all yield unphysically small values of rg. 

The results suggest that our model, with the inclusion of polarizability via 
perturbation theory, is viable for molecules with dipole moments greater than about 
3.5 D. It is of some concern that for all of the molecules studied, including those for 
which we have now found realistic values of rg, the first-order energy corrections 
are much larger than the zeroth-order energies However, one should compare the 
perturbative correction with the mean potential energy in the original dipole potential 
(Vd); not the eigenvalue (in which the negative potential energy and positive kinetic 
energy contributions noticeably cancel each other). Table 0] shows this information for 
all of the molecules studied. We see that for the most strongly polar molecules, e.g., the 
nitriles, (Vd) and are of similar magnitude. This indicates that the corresponding 
estimates of are reliable. On the other hand, for most of the other molecules, the 
magnitude of (Vd) is quite small compared to that of However, even in these 

cases the perturbation-theory estimate of the relative contribution of correlations (i.e., 
polarization) appears to be sound, at least qualitatively. 

It is interesting to compare the results from table 0] with real quantum chemistry 
calculations of positron binding to polar species. The static dipole binding energy 
is then analogous to the static, Hartree-Fock (HF) calculation of binding, while the 
total Cb can be compared with the configuration interaction (Cl) result, which includes 
correlations. In all cases the binding energy from the extensive Cl calculations is at 
least an order of magnitude greater than the HF value. For example, for hydrogen 
cyanide (HCN, p = 3D), the binding energies are 1.6meV (HF) and 35meV (Cl) [T6] : 
for formaldehyde (CH 2 O, p = 3D), 1.1 meV (HF) and 19meV (Cl) [T5|; for nitrile 
molecules (CHgCN, HCCCN, C2H3CN, C2H5CN with p = 4.1-4.4D), the HF binding 
energies are 6-18 meV, becoming 81-164 meV in the Cl calculation |19j^ . The data for 
aldehydes, ketones and nitriles in table 0] show similar large increases due to the effect 
of polarization. The model thus provides a useful estimate of the effect of correlations 
on the binding energy. 

As mentioned in the introduction, analysis of the measured binding energies found 

§ In all likelihood the above Cl energies underestimate the true binding energy, because “it is difficult 
to describe the electron-positron correlation with the established methods of computational chemistry” 

m- 
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Figure 4. Correlation between the molecular dipole moment and core radius tq, 
obtained from the perturbative polarization calculation (see Table a for nitriles 
(circles), ketones (squares), aldehydes (diamonds), formates (up triangles), acetates 
(down triangles), and methyl halides (crosses). 


the dependence of eb on a for molecules within the same chemical family (i.e., aldehydes, 
ketones, formates, acetates, nitriles) to be almost linear [7j. From (|26|) we can see that, 
for hxed p and tq, scales linearly with a. Within each chemical family, the type of 
dipole is the same and so p does not vary much. Thus, for the most part, the values of ro 
are fairly close to each other within each chemical family. This implies that considering 
the effect of polarization as the hrst-order energy correction might be quite realistic, 
even when a is large. 

4 . 2 . Dependence of the binding energy on polarizability for fixed p and tq 

The new values of tq (those obtained after including the polarizability) correlate strongly 
with the dipole moment of the molecules (see £gure|l|). However, this correlation lacks 
an obvious physical basis, and predicting the binding energy for an arbitrary molecule 
given only the values of p and a would be rather tenuous. 

On the other hand, as was stated earlier, the dipole moment p and core radius tq do 
not change vastly from molecule to molecule within each chemical family, for most of the 
families studied. To investigate the dependence of eb on a, we assigned to each family 
a hxed value of p and tq. For hve out of the six families, each with three molecules, we 
used the values of p and tq for the molecule with the median value of p. This molecule 
will hereafter be referred to as the base molecule. For the two acetates, we arbitrarily 
chose ethyl acetate as the base molecule. 

For the base molecule we know and By setting a to the appropriate values 
for the other molecules in the family, we were able to hnd the corresponding they 
are just linear rescalings of (126|) . since p and uq had not changed. With hxed by the 
values of p and tq for the base molecule, we then had estimates of Cb for every molecule 
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Table 5. The predicted values of eb found by using fixed values of y and tq for 
each chemical family, compared with the experimental values. Hydrogen cyanide is 
included with the aldehydes as it has a similar dipole moment; the value of Cb for 
hydrogen cyanide obtained using the diffusion Monte Carlo (DMC) method [T7] is also 
given. The base molecule for each family is indicated by ‘(base)’ after its name. 


Molecule 

y (D) 

ro (au) 

a (A3) 

eb 

Pred. 

(meV) 

Exp./DMC 

Aldehydes 

Butanal (base) 

2.72 

0.58 

8.2 

142 

142 

Acetaldehyde 

n 

II 

4.6 

81 

88 

Propanal 

// 

n 

6.5 

113 

118 

[Hydrogen cyanide] 

// 

n 

2.5 

45 

38 

Ketones 

Acetone (base) 

2.88 

0.63 

6.4 

172 

174 

2 -butanone 

n 

n 

8.1 

216 

194 

Cyclopentanone 

n 

n 

9.0 

240 

230 

Formates 

Propyl formate (base) 

1.89 

0.0305 

8.8 

126 

126 

Methyl formate 

n 

n 

5.1 

73 

65 

Ethyl formate 

n 

n 

6.9 

99 

103 

Acetates 

Ethyl acetate (base) 

1.78 

0.0048 

8.6 

156 

160 

Methyl acetate 

n 

n 

6.9 

125 

122 

Nitriles 

Propionitrile (base) 

4.05 

1.24 

6.3 

249 

245 

Acetonitrile 

// 

n 

4.4 

183 

180 

2 -methylpropionitrile 

n 

n 

8.1 

311 

274 

Methyl halides 

Methyl fluoride (base) 

1.86 

0.042 

2.4 

0.33 

0.3 

Methyl chloride 

// 

n 

4.4 

0.61 

25 

Methyl bromide 

// 

n 

5.6 

0.78 

40 


in the family. 

Table [5] compares the predicted and experimental valnes of the binding energy for 
the seventeen molecnles stndied. We also use this method to predict the binding energy 
of HCN, placing it in the aldehyde family (see below). The results for each family are 
also shown graphically in hgure [5l The dashed lines on the graphs show linear hts of 
the measured binding energies, while the solids lines display the linear dependence of 
the calculated binding energy on ct, as described by equation fl26|) . 

For the aldehydes, the results are very good. In particular, note the similar slopes of 
the experimental and predicted dependences of the binding energy on a. The predicted 
binding energies of acetaldehyde and propanal agree with the experimental values to 
within 8% and 5% respectively. 

For the ketones, the results are again pleasing. The predicted binding energies of 
2-butanone and cyclopentanone agree with the experimental values to within 12% and 
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Figure 5. Predicted and experimental/DMC values of eb as functions of the dipole 
polarizability. The black circles and black, solid line are the predicted values of Cb- 
The blue squares are the experimental/DMC values of Cb, with the blue, dashed line 
a linear regressive fit. 


5% respectively. 

The formates also yield good results. The predicted binding energies of methyl 
formate and ethyl formate agree with the experimental values to within 13% and 4% 
respectively. It is actually quite surprising that the predicted binding energies are as 
accurate as they are for this family, since ro for methyl formate is an order of magnitude 
smaller than the values for the rest of the family, and tq for ethyl formate is more than 
double the value for propyl formate. 

The predicted binding energy for methyl acetate is excellent; it is within 3% of the 
experimental value. Again, this is fairly surprising, given that ro for methyl acetate is 
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an order of magnitude less than the value for ethyl acetate. Note also that the a scaling 
holds well for both the formates and acetates in spite of the unphysically small values 
of Tq. 

Coming to the nitriles, we note the very encouraging results. The predicted binding 
energies of acetonitrile and 2 -methylpropionitrile agree with the experimental values to 
within 2% and 14% respectively. 

Finally, we note that for methyl halides the present model fails completely. The 
predicted binding energies for CH 3 CI and CHsBr are only 2 % of the measured values. 
At this point, we note that methyl halides are the smallest molecules examined. The 
lightest of them, methyl fluoride, also has the smallest moment of inertia, which means 
that molecular rotations neglected by the model have the largest effect on this molecule. 
This, combined with the smallest dipole polarizability, could be one of the reasons for 
the anomalously small binding energy (0.3 meV) of this molecule. Hence, when the 
“atypical” CH 3 F is chosen as the base molecule, the results for the other two molecules 
are poor. Another reason that sets methyl halides apart is that other molecules within 
each family consist of the same types of atoms. They are quite similar chemically and 
have similar ionization potentials (typically, not varying by more than 0.5 eV within 
each family [3^). On the other hand, the three methyl halide molecules contain different 
atoms (F, Cl or Br), and their ionization potentials vary considerably more: 12.47eV, 
11.22 eV and 10.54eV for CH 3 F, CH 3 CI and CH 3 Br respectively [37] . This means that 
the additional attraction due to virtual positronium formation, which is not accounted 
for by the dipole polarizability (see, e.g., m), grows along this sequence. Since this 
additional attraction is not present in our model, we obtain very poor predictions 
of binding energies. At the same time, the dipole moment, even though not much 
greater than Pcvit, plays a crucial role for binding by these molecules. Had these been 
nonpolar, atom-like species, then, based on the their ionization potentials and dipole 
polarizabilities, they would not have had bound states at all (see | 12 | for the conditions 
of binding by atoms). 

Overall, the perturbative treatment of polarization has been surprisingly good for 
the five families of organic molecules. The maximum error in any of the predicted 
binding energies is 14%, even though two of the five families exhibit small absolute 
values of tq with significant relative differences in the values of tq. The model also lends 
support to the empirical linear relationship ([ 2 ]) between the binding energy and the 
dipole polarizability, even though the values of tq for most molecules are unphysically 
small. In addition, it allows one to predict the binding energy for any molecule with a 
dipole moment that is comparable to those in one of the chemical families studied, by 
placing it in that group and rescaling using the appropriate polarizability. 

As an example, consider hydrogen cyanide (HCN), which is a linear, triatomic 
molecule with a dipole moment of 2.98 D and a polarizability of 2.5 [37]. Due to its 

toxicity, the binding energy for hydrogen cyanide has not been measured experimentally. 
Nevertheless, there already exist theoretical calculations of this energy using a variety 
of methods, such as Cl and diffusion Monte Carlo (DMC) [161 IH]- We estimated 
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the binding energy by placing hydrogen cyanide among the aldehydes, since they have 
similar dipole moments. The resulting prediction of the binding energy (45meV) is 
within 18% of the DMC value of 38meV [H] (see table [5] and hgure [5]). Note that 
the ketones actually have more similar dipole moments to hydrogen cyanide than the 
aldehydes. Placing HCN in the ketone family led to a predicted binding energy of 
70 meV, which is a factor of two greater than the DMC result of HU. This large value 
is likely an overestimate, in spite of the fact that quantum chemistry calculations tend 
to give lower bounds for the positron binding energies [ini m\. Experimental data 
for ketones shows signihcant deviations from linearity (see hgure [5]) , which makes the 
ketone-based prediction for HCN less relaible. 


5. Non-perturbative treatment of molecular polarization 


5.1. Models of polarization potential and core radii 


Although the perturbative inclusion of polarization described in section |4] generates 
reasonably accurate predictions of binding energies for organic molecules, the effective 
core radii tq are too small for most molecules to be physically meaningful. In addition, 
the hrst-order polarization energy correction for most molecules is too large to justify 
the use of perturbation theory. To overcome this limitation, in this section we include 
the polarization potential in a full, non-perturbative manner. As we will see, this leads 
to new physical insights and hnally gives good physical values of tq for all molecules. 

When the polarization potential fl2T|) is added to the Schrodinger equation (0]), the 
angular equation (|^ remains unchanged, and we have the radial equation 


1 d^P 


2 dr^ 


+ 


2r2 


2r‘^ 


Pmsnij'^ EmsnPmsni. 


(27) 


for the function Pmsuir) = f'Rmsnij')- Here we have also introduced the polarization 
cut-off function g{r), which tends to unity at large r and moderates the unbounded, 
unphysical growth of the —a/2r^ term at small distances (see, e.g., [3]; also see below). 

At large r, the polarization potential is negligible in comparison with the 1/r^ dipole 
potential. Thus, at some sufficiently large value of r = r ma ^ the radial wave function is 
given by equation (fTT|) . which gives the boundary conditions 


P 


msny max 

d-Pmsnij ') 


dr 


^msn \/^maxA"i/3ms max) •> 

^msn \p/^ K\l3rns^^''TLSrJ' 


(28 a) 
(286) 


where Amsn is an arbitrary constant. A value of r^ax = 30 an has been used throughout. 
Solving equation fl27)) numerically in the interval 0 < r < r^ax with m = 0, s = n = 1, 
Poll = —Cb, and Aon = ''"max yields a real function Poii(r) with infinitely many roots 
accumulating at r = 0. The largest of these roots is the value of tq. 

We initially considered g{r) = 1, as in section 01 This led to values of tq in the 
range 1.55-2.32 an across the six families of molecules, which are much greater than 
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their perturbative counterparts, and probably too large to be considered physical. This 
is due to the polarization potential fl2T|) blowing up and causing a rapid variation of the 
radial wave function at small r, while in reality the polarization is a long-range effect. 
For the same reason, the binding energy was found to be extremely sensitive to the 
value of ro, making it very difficult to use the model in a predictive way. 

Consequently, we considered several cut-off functions, viz.. 


9i{r) 

= l-exp(-rVr6), 

(29 a) 


^4 

(29&) 

92{r) 

(r+ rc)4’ 

93{r) 


(29 c) 

(r2 -K r2)2’ 


where Uc is a cut-off radius for the polarization. The function gi{r) provides a very 
rapid cut-off, and is commonly used to model polarization potentials in atoms [3]. The 
functions and g^{r) vary much more slowly. They can effectively account for the 
fact that the “centre” of the polarization potential is usually off-set with respect to the 
location of the molecular dipole. The dipole moment is usually associated with one of 
the terminal bonds, which is near one of the “ends” of the molecule rather than in the 
middle. Initially we worked with fixed values of Tc across the entire set of molecules, 
and though this reduced the values of tq from those of the “hard” potential {g{r) = 1), 
and also reduced the sensitivity of the binding energy to ro, it did not signihcantly 
reduce the large spread in tq within or between families. This led us to consider using 
a polarizability-dependent cut-off radius. 

The polarizability of organic molecules is generally proportional to the number of 
atoms or number of bonds in the molecule. This idea is the physical basis behind various 
additivity methods for the calculation of molecular polarizabilities [13]. In this spirit, 
the polarization potential at large distances is the sum of terms —aj/2r^ due to the 
contribution of individual atoms or bonds z, with the distance rj measured accordingly. 
In the spherically-averaged form fl2T|) . the distance r must measured from the “centre 
of polarization” rather than the centre of the molecular dipole. As a result, at small r 
the singular form —a/2r^ must be replaced by a constant —a/2rf as described by the 
cut-off functions ( 72 and Qsir). Here is the effective radius of the molecule. It is 
physical to link it to the polarizability a by, e.g., 

Tc = Ca^, (30) 

with C being an adjustable parameter. The polarizability a has dimensions of volume 
(i.e., length cubed), so the choice p = ^ would be sensible for three-dimensional, 
approximately spherical molecules, while /i = | would be better for approximately 
planar molecules. Experimentation showed that p = ^ works best for our set of 
molecules, with C being chosen separately for each family to minimize the range of 
values of tq within the family. 

Table |6] shows the final values of uq obtained for the molecules, with As 

expected, the cut-off functions g^if) and gsi^r) gave the most physically meaningful 
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Table 6. Values of rg when a soft polarization potential with the cut-off function 
53 (r) is included non-perturbatively and Tc given by equation (1501) with /i = i. The 
parameter C is given in units of where oq is the Bohr radius. 


Molecule 

m(d) 

a (A^) Cb (meV) 

To (au) 

Aldehydes (C = 1.08) 





Acetaldehyde 

2.75 

4.6 

88 

1.17 

Propanal 

2.52 

6.5 

118 

1.09 

But anal 

2.72 

8.2 

142 

1.16 

Ketones (C = 0.98) 





Acetone 

2.88 

6.4 

174 

1.24 

2 -butanone 

2.78 

8.1 

194 

1.24 

Cyclopentanone 

3.30 

9.0 

230 

1.37 

Formates (C = 1.06) 





Methyl formate 

1.77 

5.1 

65 

0.94 

Ethyl formate 

1.98 

6.9 

103 

0.98 

Propyl formate 

1.89 

8.8 

126 

0.94 

Acetates (C = 0.96) 





Methyl acetate 

1.72 

6.9 

122 

1.05 

Ethyl acetate 

1.78 

8.6 

160 

1.05 

Nitriles (C = 1.12) 





Acetonitrile 

3.93 

4.4 

180 

1.34 

Propionitrile 

4.05 

6.3 

245 

1.34 

2 -methylpropionitrile 

4.29 

8.1 

274 

1.43 

Methyl halides (C = 0.95) 





Methyl fluoride 

1.86 

2.4 

0.3 

1.22 

Methyl chloride 

1.90 

4.4 

25 

1.24 

Methyl bromide 

1.82 

5.6 

40 

1.24 


results, with neither being signihcantly better or worse than the other. Here we 
present the results obtained using g 3 {r). Figure [6] compares the radial function Poii(^) 
for acetonitrile, with and without the non-perturbative inclusion of the polarization 
potential. It is clear from the hgure that the wave function does not change much for 
r > 5 an. However, at smaller distances the addition of the polarization potential causes 
a more rapid variation of the wave function, leading to a greater core radius tq. 

As seen from table |6l all of the radii are now ~1 au and hence look physically 
meaningful. They also remain approximately constant within each chemical family. 
The maximum range of tq within any family is 0.13 au (for the ketones), and the range 
across all the molecules is 0.49 au. The values of C are also quite consistent, ranging 
from 0.95 to 1.12 oq As before, the largest radii are obtained for the most strongly 

polar molecules (the nitriles) and the smallest radii are obtained for the most weakly 
polar molecules (the formates and acetates). 
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Figure 6. Radial wave functions for acetonitrile. The dashed black curve is without 
the inclusion of polarization; the solid blue curve is with the polarization included 
non-perturbatively, using the cut-off function 33 (r). 


5.2. Dependence of the binding energy on polarizability for fixed p and tq 

We can now again investigate the dependence of the binding energy on the molecular 
polarizability, by hxing p and rg within each family and varying a. In these calculations 
we choose the same base molecule within each family as in section 14.21 The binding 
energy which enters in equation fl27|) is then adjusted for each molecule in the family 
until the core radius rg of the base molecule is obtained. Note that the polarization 
potential is now included non-perturbatively, hence, there is no reason to expect that eb 
depends lineary on a. The resulting binding energies are shown in table [7] and hgureO 

There is generally very close agreement between the model predictions and the 
experimental data; for every molecule except methyl formate and methyl acetate, 
the relative difference of the predicted binding energy from the experimental value 
is smaller than it was using the perturbative method. Particularly noteworthy is 2- 
methylpropionitrile, for which the predicted binding energy coincides exactly with the 
experimental value. The error for methyl formate has increased from 13% to 18%, and 
for methyl acetate it has increased from 3% to 4%. 

From hgure [7] it is apparent that when the polarization potential is included non- 
perturbatively, the dependence of eg on a is indeed non-linear. For all families except 
the methyl halides, Cb increases convexly with a for a < 4A^; for a > 4A^ the growth 
becomes concave. The molecules studied all lie in the convcave region, and the growth 
for them could be reasonably well approximated by a straight line. The prediction 
curve for the methyl halides is different but particularly remarkable, as the description 
of this molecular family was extremely poor in the perturbative treatment. Here the 
molecules lie in the convex region (which spans a larger range of polarizabilities than 
for the other families), and the dependence of Cb on a is markedly nonlinear. However, 
close agreement is observed with the measured binding energies for CH 3 CI and CHgBr. 
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Table 7. The predicted values of eb found by using fixed values of y and tq for each 
chemical family with the polarization included non-perturbatively, compared with the 
experimental values. Hydrogen cyanide is included with the aldehydes as it has a 
similar dipole moment; the value of eb for hydrogen cyanide obtained using the diffusion 
Monte Carlo method m is also given. The base molecule for each family is indicated 
by ‘(base)’ after its name. 


Molecule 

y (D) 

ro (au) 

a (A3) 

eb 

Pred. 

(meV) 

Exp./DMC 

Aldehydes 

Butanal (base) 

2.72 

1.16 

8.2 

142 

142 

Acetaldehyde 

// 

n 

4.6 

86 

88 

Propanal 

// 

ff 

6.5 

122 

118 

[Hydrogen cyanide] 

n 

tf 

2.5 

32 

38 

Ketones 

Acetone (base) 

2.88 

1.24 

6.4 

174 

174 

2 -butanone 

n 

tf 

8.1 

210 

194 

Cyclopentanone 

n 

tf 

9.0 

225 

230 

Formates 

Propyl formate (base) 

1.89 

0.94 

8.8 

126 

126 

Methyl formate 

// 

tf 

5.1 

77 

65 

Ethyl formate 

// 

tf 

6.9 

106 

103 

Acetates 

Ethyl acetate (base) 

1.78 

1.05 

8.6 

160 

160 

Methyl acetate 

// 

tf 

6.9 

127 

122 

Nitriles 

Propionitrile (base) 

4.05 

1.34 

6.3 

245 

245 

Acetonitrile 

// 

ff 

4.4 

200 

180 

2 -methylpropionitrile 

// 

ff 

8.1 

274 

274 

Methyl halides 

Methyl fluoride (base) 

1.86 

1.24 

2.4 

0.3 

0.3 

Methyl chloride 

// 

ff 

4.4 

20 

25 

Methyl bromide 

// 

ff 

5.6 

42 

40 


Table [7] and figure [7] also show our estimate of the binding energy for HCN. As 
seen from the graph, the experimental value lies very close to the prediction curve for 
the aldehydes data, which provides support that the dependence of Cb on a is not truly 
linear. The value obtained (32meV) is within 16% of the DMC calculation (38meV) 
na, which is slightly closer than our perturbative estimate. Interestingly, if we now 
place HCN in the ketones family, its estimated binding energy becomes 40 meV, in 
very close accord with the DMC value. This is further evidence that a nonperturbative 
treatment of molecular polarization gives overall much more consistent results. 
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Figure 7. Predicted and experimental/DMC values of es as functions of the 
dipole polarizability, obtained using the non-perturbative inclusion of the polarization 
potentil. The black circles are the predicted values of Ch, with the black, solid curve 
showing the calculated dependence of Cb on a for each family. The blue squares are 
the experimental/DMC values of Cb, with the blue, dashed line a linear regressive fit. 


6. Concluding remarks 

Here we provided a simple model for positron binding to polar molecules, which captures 
the essential physics of this system. 

Modelling the molecule as a sphere of radius ro with a static point dipole of dipole 
moment p at the centre and using experimental data on binding energies required 
unphysically small values of tq, even for the most strongly polar molecules. This 
indicated that the binding energies are greatly enhanced by some factor other than 
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the molecule’s permanent dipole moment, i.e., electron-positron correlations. 

Including the effect of correlations perturbatively through the polarization potential 
did confirm this expectation. It showed that even for the strongly polar molecules, the 
effect of correlations increased the binding energy by an order of magnitude compared to 
the static-dipole calculation. The observed increase matched the difference between the 
Cl and HF binding energies obtained in state-of-the-art quantum chemistry calculations. 

Including the polarization potential as a perturbation of the original Hamiltonian 
also yielded larger, more physical values of ro for all of the molecules studied, but for 
most molecules they were still too small to be interpreted directly. This was partly 
due to the fact that the true static potential for the positron near the molecule is less 
repulsive than the hard wall of our model. Reduced values of tq may also account 
for some of the short-range correlation effects, such as virtual positronium formation. 
Sensible values of ro were, however, obtained for the nitriles, the most strongly polar of 
all the molecules studied. In spite of the fact that most of the molecules had unphysical 
values of ro, it was found that taking the value of p and ro for the molecule in each 
chemical family with the median dipole moment and varying the polarizability to match 
the other molecules in the family, gave reliable predictions of the binding energies for 
those molecules (with the exception of the methyl halides). The perturbative treatment 
was also in line with the observation made by experimentalists, that the dependence 
of the binding energy on the polarizability of the molecule is apparently almost linear 
[7]. Of course, a general increase in the binding energy with the polarizability could 
be expected, but there was no explanation for the linear dependence. According to 
our model, this feature indicates that the perturbation theory is at least qualitatively 
correct, even though the first-order energy corrections are generally greater than the 
original (dipole) eigenenergies. The results of this treatment for the methyl halides, 
however, were very poor in comparison with the other families, and we attributed this 
to the fact that the binding by the base molecule (CH 3 F) is likely affected by rotations, 
and that the three methyl halide molecules had a signihcantly larger range of ionization 
potentials that the other families. The latter could indicate a significant change in the 
contribution of virtual positronium formation across the members of this family, which 
is not accounted for in our model. 

A test of the model was to use it to predict the binding energy for hydrogen cyanide, 
which has not been measured experimentally. Our estimate agreed with a previous 
calculation using the diffusion Monte Carlo method to within 20%, which provides 
evidence that our model has good predictive power and could be useful for estimating 
the binding energies that have never been measured in experiment, provided that the 
binding energy for a molecule with a similar value of p is known. 

The most glaring limitation of the model as it stood was that it could not predict 
binding energies using only the dipole moment and polarizability of a molecule. To 
perform a calculation one needed a value of tq, and unless binding energies for molecules 
with similar values of p were known, one could not easily choose a suitable value for 
tq. In fact, we found that for most of the chemical families we considered the values 
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of ro had no immediate physical relevance, and for the most weakly polar families (i.e., 
formates and acetates) there were signihcant variations in the values of tq despite the 
similar values of p. 

In a bid to attain physically meaningful values of tq for all of the molecules, 
we proceeded to include the effects of polarization in a non-perturbative way. We 
experimented by solving the radial Schrodinger equation numerically using several 
model polarization potentials, and found that the best results were obtained using a 
polarization potential of the form —a/2(r^ + r^)^, with the cut-off radius Vc = 

[C being a constant). The parameter C was chosen separately for each chemical 
family so as to minimize the spread of ro within each family. Physically meaningful 
values of tq (in the range 0.94-1.43 au) were then obtained for all molecules. By 
choosing the values of C carefully, the spread of values of ro within in each family 
was made relatively small, though there was inevitably still a larger range of 0.49 au 
across the entire set of molecules. The most strongly (weakly) polar molecules still 
possessed the largest (smallest) values of ro. Again hxing p and ro for each family 
and varying only the polarizability led to excellent predictions of Cb; the predictions 
were generally more accurate than their perturbative counterparts, with particularly 
large improvement for the methyl halides. We observed that the true dependence of 
Cb on a is actually nonlinear. The prediction for HCN was also slightly better than its 
perturbative counterpart, and supported our observation of nonlinear growth of Cb with 
a. 

In summary, our model can be used effectively to predict positron-molecule binding 
energies based on the molecular dipole moment and dipole polarizability, particularly 
when polarization is included in a non-perturbative way. It provides a clear picture of 
the system, thereby complementing the current computational effort towards rigorous 
theory of positron-molecule binding. 
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